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ABSTRACT 

The Sun is believed to have been the recipient of a substantial amount of 
metal-rich material over the course of its evolution, particularly in the early 
stages of the Solar System. With a long diffusion timescale, the majority of 
this accreted matter should still exist in the solar convection zone, enhancing 
its observed surface abundance, and implying a lower- abundance core. While 
helioseismology rules out solar models with near-zero metallicity cores, some 
solar models with enhanced metallicity in the convection zone might be viable, 
as small perturbations to the standard model. Because of the reduced interior 
opacity and core temperature, the neutrino flux predicted for such models is 
lower than that predicted by the standard solar model. This paper examines 
how compatible inhomogeneous solar models of this kind are with the observed 
low and intermediate degree p-mode oscillation data, and with the solar neutrino 
data from the SNO Collaboration. We set an upper limit on how much metal-rich 
accretion took place during the early evolution of the Sun at ~ 2M^ of iron (or 
~ 40M^ of meteoric material) . 

Subject headings: Sun: accretion, evolution, interior, oscillations, neutrinos 



1. Introduction 



Early on in the evolution of the Sun, the interstellar medium and young Solar System 
provided an environment capable of frequent accretion events. Material in the form of ISM, 
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proto-planets, planetesimals, comets, and asteroids frequently bombarded the Solar surface 
and were accreted. This accreted matter is deficient in hydrogen and helium, introducing 
metal-rich material into the upper layers of the Sun. Exactly how much metal-rich matter 
was accreted in these early stages and for how long, has been speculated by several authors 
(Christenscn-Dalsgaard, Gough, & Morgan 1979, Jeffery et al. 1997), but unfortunately is 
not well known. 

More recently, Murray et al. (2001) have searched a sample of 640 solar-type stars for 
the signature of iron enhancement in their spectra and concluded that, on average, these 
stars appear to have accreted about O.5M0 of iron while on the main sequence. They raise 
the possibility that the Sun may have accreted a similar amount of iron during its evolution. 
Whether helioseismology can detect this amount of accretion is uncertain. So far only one 
study has been published on the effect of accretion on the p-mode frequency (Henney & 
Ulrich 1998). These authors concluded that the accretion of SM® of meteoric material 
(approximately O.4M0 of iron) on the Sun could not be detected by seismology due to other 
uncertainties in the models. 

The main consequence of such accretion is the metal enrichment of the Sun's surface and 
convection zone, while maintaining a lower-abundance core. In such a scenario, the expected 
neutrino flux is reduced due to the lower opacity and temperature in the central region of the 
Sun. It has long been known that solar models with near-zero metallicity cores can lower the 
predicted neutrino flux, and this explanation was discussed early-on as a possible sohition of 
the classical neutrino problem (sec e.g. the review by Rood 1978). More recently, Gucnther 
& Demarque (1997, hereafter GD97) have constructed solar models with low-Z cores using 
present physics and found their p-mode frequencies to be incompatible with solar p-mode 
observations. 

This paper describes work which can be viewed as an extension of GD97; but this 
time we examine the limits set by helioseismology on the non-standard assumption that the 
solar convection zone of the Sun has been enriched in heavy elements by accretion during 
its early evolution. In addition to comparing the p-mode frequencies to observation in the 
frequency difference diagram, as done by GD97, we also compare the calculated sound-speed 
and density in our models to the same quantities derived from observation by inversion (Basu 
& Antia 1997, hereafter BA97; Basu, Pinsonneault & Bahcall 2000, hereafter BPBOO). 

Since this research was begun, the first results from the Sudbury Neutrino Observatory 
(SNO) experiment have been released (Ahmad et al. 2001). Ahmad et al. (2001) conclude 
that the ^B neutrino flux derived from their observations, i.e. $8^ = 5.44 ± 0.99 x 10^ 

cm~^s~^, agrees with the predictions of the best standard solar models. Any accretion model 
must then also be constrained by the SNO observations, in addition to the constraints of 
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helioseismology. 

2. Assumptions and Method 

The main assumption in this paper is that the accretion of heavy elements onto the 
Sun took place during the early evolution of the Sun, at a time when the Sun was near the 
main sequence, i.e. when the convection zone had reached its main sequence value. The 
convection is shallow at this time, and only the outer 2% of the Sun's mass in the convection 
zone get enriched. 

As in GD97, the solar models were constructed using the Yale stellar evolution code 
(YREC for Yale Rotating Evolution Code) in its non-rotating configuration (Guenther et al. 
1992). The energy generation routines used in YREC are those of Bahcall & Pinsonneault 

(1992). The neutrino cross sections used are the same as in Bahcall, Pinsonneault & Basu 
(2001). The equation of state tables prepared by the OPAL researchers (Rogers 1986, Rogers, 
Swenson, & Iglesias 1996) were used for the models. The interior opacities were derived 
from the OPAL opacity tables (Iglesias & Rogers 1996), while the surface and atmospheric 
opacities were taken from the Alexander & Ferguson (1994) tables. 

All models were evolved from a zero-age main sequence (ZAMS) model to near the age 
of the Sun in 50 equally spaced time steps. The age of the Sun was taken to be 4.5 Gyr, 
close to the meteoric solar age (Guenther 1989). As was noted by GD97, standard solar 
models of this age are in best agreement with the calculated oscillation spectra. The mixing 
length/pressure scale height ratio, a, and the helium content by mass, Y, were adjusted 
automatically by the code to produce a precisely calibrated standard solar model (Guenther 
et al. 1992). Most models were tuned to match the observed solar radius to one part in 
10^ (R0=6.958 X 10^°cm) and the solar luminosity to one part in 10^ (L0=3.8515 x lO^^ergs 
s~^). Models #15, 16, and 19 (see Table 1) did not converge at the aforementioned tolerance 
for radius, but converged to the solar radius to one part in 10®. The number of shells in 
standard solar models #17 and 18 was 1892, and 1893 in models #19 and 20. There were 
more than 2300 shells in all non-standard solar models. Tests have shown that at least that 
many shells are needed to obtain the required precision in the calculation of the p-mode 
frequencies (Guenther et al. 1992; Guenther 1994). 

Helium and heavy element diffusion were included in all of the solar models due to their 
necessary role in finding simultaneous agreement with both the observed value of (Z/X)© 
and the p-mode spectrum (for further details see GD97). The assumption made regarding 
the nonstandard models is that due to the accretion of metal-rich matter, the Sun's interior. 
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(defined here as R < Renv), is composed of material at a lower heavy-element abundance than 
the Sun's surface. In order to produce such nonstandard solar models with low-Z interiors, 
the run of Z in the initial ZAMS model was modified. The interior metal abundance was 
initially set to the homogeneous value of Zi^t out to (M/M0)=O.9, indicating a metal-poor 
interior. The more metal-rich exterior of the Sun, (M/M0)>O.975, was setup with a metal 
abundance of Zinit- Zint and Zinit represent ZAMS mass fractions of all heavy elements for 
the interior and exterior, respectively. In the intermediate region, O.9<(M/M0)<O.975, Z 
linearly increases from Zint to Zinit. 

In deciding what value Zint might be for the nonstandard solar models, it is relevant 
to note that the Sun is observed to be more metal rich than the surrounding ISM, with 
ZisM possibly as low as 65% of Zq (Mathis 1996). Beheving the Sun formed from material 
typical of the ISM, the metal-enhanced exterior could have resulted from the bombardment 
of metal-rich material in the form of comets, asteroids, planetesimals and proto-planets. 
With a long diffusion timcscale, most of this material should still exist in the upper layers 
of the Sun, leaving the interior metal-poor, much like the surrounding ISM. With this in 
mind, we chose a value of Zint=0.65ZQ. Other values of 0.30Zq, O.SOZ©, and O.SOZ© were 
also examined. 

A choice of Zint = O.SOZinit corresponds to an accretion enhancement of about 2M0 
in iron, or about 4OM0 in meteoric material. Similarly, Zint=0.65Zinit, Zint=0.50Zinit, and 
Zint=0.30Zinit corresponds to about 2.9M0, 4.IM0 and 5.7M0, respectively, in iron accretion, 
or about 6OM0, SOM® and I2OM0, respectively, in meteoric material accretion. 

Four standard solar models were created, differing only by the value of Zinit, the initial 
or ZAMS mass fraction value of all heavy elements in the solar exterior. Without the 
assumption of a more metal-poor interior, Z^nt is assumed equivalent to Zinit. Values of Zinit 
were taken as 0.0170, 0.0188, 0.0200, and 0.0220. These standard solar models comprise 
models #17-20 in Tables 1-3. Zinit was also varied in the nonstandard solar models, thus 
along with the varying Zint, a grid was created, comprising models #1-16 in Tables 1-3. 

Physical characteristics of both the standard (#17-20) and nonstandard (#1-16) solar 
models are listed in Table 1. Table 1 includes, from left to right: Model, the model number; 
Type, the type of model, where NSSM stands for a nonstandard solar model and SSM stands 
for a standard solar model; Xinit, the initial or ZAMS mass fraction value of hydrogen; Zinit, 
the initial or ZAMS mass fraction value of all heavy elements in the solar exterior; Zint /Zinit, 
the initial or ZAMS mass fraction value of all heavy elements in the solar interior relative 
to Zinit; Xgurf, the surface mass fraction value of hydrogen at the evolved age; Zgurf) the 
surface mass fraction of all heavy elements at the evolved age; Mcnv, the fraction of the total 
mass contained in the outer convective envelope; Renv, the radius fraction of the base of 
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the convective envelope; log Pc, the base ten logarithm of the central pressure; logTc, the 
base ten logarithm of the central temperature; log p^-, the base ten logarithm of the central 
density; Xc, the central mass fraction of heavy elements; Zc, the central mass fraction of 
heavy elements. 

The nuclear energy generation properties of both the standard (#17-20) and nonstan- 
dard (#1-16) solar models are listed in Tabic 2. Tabic 2 lists the fraction of total photon 
luminosity coming from the PP 1, PP II, and PP III branches of the PP network and from 
the CNO cycle. Also listed are the individual neutrino fluxes from the neutrino producing 
reactions that occur in the sun (see Figure 1 of GD97). Note that in GD97 the neutrino 
fluxes are listed in units of 10^°cm~^s~^. Finally, Table 2 lists $(^^01), the total neutrino 
flux, in SNU, for the ^'''Cl detector; and $(^^Ga), the total neutrino flux, in SNU, for the 
^^Ga detector. 



3. p-mode Frequencies 

As a test on the validity of a model, the expected oscillation frequencies can be compared 
with observations from the Sun. Guenther's non-radial, non-adiabatic pulsation program 
(Guenther 1994) was used to calculate the oscillation frequencies of the models produced 
with YREC. The model output was then compared with data obtained by the Michelson 
Doppler Imager (MDI) instrument on board the Solar and Heliospheric Observatory (SOHO) 
during the flrst year of its operation (Rhodes et al. 1997). This data set was chosen as it 
comprises one of the longest time-series of data, 360 days. More recent MDI-SOHO data 
include only 144-day or 72-day data sets, and therefore restrict the number of data points, 
with most sets having very few /=0,1,2 modes. Data from the GOLF experiment (see Bertello 
et al. 2000 and Garcia et al. 2001 for the latest results) were not considered here despite the 
low-degree low-frequency modes. While these modes can tell us a lot about the structure 
of the core (see Turck-Chieze et al. 2001), for our purposes here, it requires mixing one 
set of low degree modes with intermediate and high degree modes from another instrument, 
which can give rise to artifacts in the solar core unless the data are contemporaneous (see 
Basu et al. 1996, 1997). Hence, in an attempt to avoid such artifacts, we have opted 
for a homogeneous set of data. For comparison, BiSON+LOWL (Basu et al. 1997) data 
were also used. This data set had been specially constructed by obtaining frequencies from 
contemporaneous BiSON and LOWL time scries. This set gave results very similar to those 
obtained by the MDI-SOHO data set used here. 

From these data, p-moAe frequency differences: fmodci-'^oj were computed for /=0-4, 10, 
20, 30, 40, 50, 60, 70, 80, 90, and 100. Frequency difference plots (model frequency minus 
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observed frequency versus observed frequency) were then constructed for each model and 
can be seen in Figure 1. Each hne joins together data of a common /-value, hence joining 
together j»-modes with approximately similar inner turning points. In the event of perfect 
seismic agreement with the Sun, a frequency difference plot would show a straight horizontal 
line at fiHz, indicating that the p-mode frequency differences were zero. Instead, Figure 1 
shows a more complex bundle of lines, indicating discrepancies of our models from the Sun. 
As in GD97, the quality of a models' agreement with the Sun is judged on the tightness of 
the bundle of lines. The underlying slope error present is thought to be due to modeling 
errors in the very outermost layers of the Sun, a region of known uncertainties. The errors 
in the interior of the model, primarily near the base of the convection zone, are directly 
correlated to the bundle thickness in the sense that the tighter the bundle, the better the 
model is at fitting the region near and above the convection zone. (See GD97 for further 
details.) 

Figure 1 is a 4x5 grid of plots showing the frequency differences for each of the 20 
models computed. Each plot is annotated by its model number (see Tables 1-3), the surface 
Z/X ratio, and neutrino flux, in cm~^s~^. Lines connect common /-values, with 1=0-4, 
10, 20, 30, 40, 50, 60, 70, 80, 90, and 100. As indicated by the model numbers, the four 
standard solar models computed in this work comprise the last row in Figure 1. Figure 2 
shows a larger version of the frequency-difference plot for model #20, including error bars 
from the MDI-SOHO data. The error bars indicate lOcr errors in the data averaged over 
250/xHz frequency bins, for all /-values included in this work. As can be seen in Figure 
2, the errors increase with increasing frequency. Despite the relatively good agreement of 
our model with the Sun, the accuracy of the observations is such that our best model is 
discrepant at roughly the 30cr level at low frequencies, and is worse at higher frequencies. 
This large discrepancy seen in the calculated frequencies is not found in the sound speed 
comparison (sec Figure 6). It is primarily due to modeling errors in the surface layers of the 
Sun, which could be reduced by further analyses including magnetic fields, turbulence, and 
a better understanding of the convection zone. 

Special consideration was given to low / j)-modes due to their ability to probe the deep 
interior of the Sun. This is due to the fact that the inner turning points of these modes are 
located closest to the core. Low / p-modes are still sensitive to the outer layers, but this effect 
can be canceled out by subtracting from a given |?-mode frequency the frequency of a |?-mode 
with a similar eigenfunction shape in the outer layers and distinct eigenfunction shape in the 
deeper layers. The small spacing difference, defined as 5u{n, I) = ^{n, I) — v{n — 1, /+2), thus 
provides a further diagnostic of the deep interior of the Sun. (See GD97 for further details.) 
To compare our models with observations from the Sun, small spacing differences plots were 
produced [i.e., SumoAei^n, l) — 5uQ{n, I) vs. ^'©(n, /)]. As in the case of the frequency difference 
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plots, the closer 5^'modei('^5 ~ ^^©('^j is /uHz, the better the agreement between the 
observations and a model. 

Since the sensitivity of the small spacings to the deep interior of the Sun diminishes with 
increasing I, only those p-modes with /=0, 1, and 2 were used to contrast the models with 
observed values from the Sun. The small spacing difference plots can be seen in Figure 3. 
Figure 3 is a 3x5 grid in which each row contains models of the same Zint, where Zint=0.30Zinit 
includes models #1-4, Ziiit=0.50Ziiiit includes models #5-8, Zint=0.65Zinit includes models 
#9-12, Zint=0.80Zinit includes models #13-16, and Zint=Zinit includes the standard solar 
models #17-20. Each column in Figure 3 contains results for a common /-value, with /=0 in 
column 1, 1=1 in column 2, and 1=2 in column 3. Lines in each panel represent models with 
the same Zjnt with differing Zi^it values ranging from 0.0170 to 0.0220, as indicated. 

4. Sound-speed and Density 

Figure 4 shows the relative sound-speed differences for the twenty models presented in 
this paper compared to the observed solar sound-speed, derived from an inversion of the solar 
p-mode frequencies (BPBOO). We note the superior agreement of the standard solar models 
with observation, especially model #20, as compared to the non-standard models. Similarly, 
Figure 5 shows the same plots for density. The less precise agreement with densities than 
with sound-speed is in part a reflection of the greater uncertainty in density inversions. 

5. Solar Models 

5.1. Choice of a Benchmark Standard Solar Model 

Figure 1 shows the characteristics of the standard solar models (#17-20) in the frequency 
difference diagram. Best agreement with observation, as measured by the thickness of the 
line bundles corresponding to different /-values, favors #20 and 19. Note that our model #19 
is nearly identical to model #20 of GD97, that was judged by GD97 as the "best" standard 
solar model in their study. But our model #20, which is intermediate in Zinit between GD97's 
models #20 and 21, is a better model still. The bundle of /-values in model #20 is thinner 
than in model #19, a fact which is also reflected in the best agreement with observation in 
the sound-speed plot shown in Figure 4. 

In selecting a "best" standard solar model, one should also take into account two addi- 
tional constraints. The first constraint, to which we will assign the most weight, is the solar 
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convection zone depth, Renv, which has been derived by sound-speed inversion (Christensen- 
Dalsgaard et al. 1991; BA97). The second additional constraint is Ysurf for the Sun. This 
quantity has also been estimated by inversion of the p-mode data (Basu 1998). 

For model #20, we have Renv = 0.7128, to be compared to Renv = 0.713 ± 0.001 from 
inversion (BA97). This agreement is confirmed by the sound-speed plot for model #20 (see 
Figures 4 and 7), which shows close agreement with the inversion data. At the same time, 
model #20 yields Ygurf = 0.2519, to be compared to the helium mass fraction 0.248 ± 0.003 
derived by Basu (1998), a satisfactory agreement in view of the uncertainties in diffusion 
coefficients and efficiency. In addition, the ^B neutrino flux of 5.63 cm~^s~^ for model #20 
is in excellent agreement with $8^ = 5.44 ± 0.99 x 10^ cm~^s~^ found by SNO (Ahmad 
et al. 2001). On the other hand, the surface metaUicity (Z/X)surf is a little too large. One 
can see that a model interpolated between models #20 and 19 would satisfy all constraints 
within the errors. 

As it is, model #20 is an excellent standard solar model, as illustrated in Figures 4-8. 
Comparing Figures 2 and 3 in BPBOO with our Figures 7 and 8, we see that model #20 
compares favorably with model BP2000. 

The small spacing plots (Figure 3) shows that all our standard models, models #17- 
20, agree equally well with observation. Within the range of chemical composition we have 
considered, small spacings are insensitive to the choice of initial heavy element content. 
Figure 6 shows the zoomed-in small spacing plot for l—O for standard solar model #20. 
Error bars calculated from the MDI-SOHO data set used are also plotted. The observational 
uncertainties are approximately ±0.07 //Hz for the small spacing determination. From this 
figure, the numerical uncertainties in the solar model appear to be approximately ±0.05 to 
0.10 iiHz for the small spacing calculation. The larger uncertainties of ±1 fiHz for the mode 
frequencies themselves cancel out somewhat when calculating the small spacings. 

It is difficult to evaluate the uncertainty in the observed solar heavy element content 
(Z/X)surf. GD97 used (Z/X)surf = 0.0244 ± 0.001, from Grevesse, Noels & Sauval (1996). 
More recent lower estimates suggest that large systematic errors may still exist in the mea- 
surements. Grevesse & Sauval (1998) find (Z/X)surf = 0.0230, quoting a 10% error estimate. 
Asplund (2000), using a stellar model atmosphere constructed with a more realistic treat- 
ment of convection, derives (Z/X)surf = 0.0226. Unfortunately, theoretical estimates for both 
Ygurf and (Z/X)su].f are also difficult to make, due to the uncertainties in calculating the ef- 
ficiency of diffusion just below the convection zone in the stellar interior (Chaboyer et al. 
1995). Helioscismic sound-speed inversion shows that the tachocline structure is one of the 
least well-understood features of standard solar models (BA97). 
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5.2. Non- Standard Solar Models 

At first glance, the frequency difference diagram in Figure f does not reveal much 
difference in bundle thidcness between the non-standard models ^^7, 8, ff, and f5 and the 
standard model #19. The /-value bundle for model #12 seems even slightly thinner. But 
models #7, 8, 11 and 12 are all readily ruled out by the neutrino constraint. 

Non-standard model #16 is similar to standard model #20 in having the thinnest l- 
value bundles and a marginal value for (Z/X)surf- But it fails the ^B neutrino test by a small 
amount. Non-standard models #10 and 15, like standard model #20, pass the convection 
zone depth test, but of these, only #15 passes the (Z/X)surf test. However, model #15 fails 
the ^B neutrino test. 

Turning to the small spacing diagram (Figure 3), we see that the non-standard models 
agree as well with observation as the standard models only for / = 0. But for / = 1 and I — 
2, there is an increasing discrepancy with decreasing Zint. The maximum discrepancy is only 
2 fjHz, however, which may still be within the uncertainties. As in the case of the standard 
models, heavy element content plays little role. 

As a further comparison of the models calculated here, the root-mean-square of the 
sound-speed difference, rms(5c/c), was computed for each model. Seen in Table 3, the results 
indicate that those models with Zinit=0.0220 fair well. Model #20 certainly outperforms the 
others, followed by standard solar model #19 and non-standard solar model #16. These 
results further strengthen the argument that model #20 is in best agreement. Table 3 
summarizes the values of the constraints considered here for each solar model computed in 
this work. Listed are the Model, Renv) Ygurf, (Z/X)surf and rms((5c/c). 

We conclude that of all the non-standard solar models listed in Tables 1 and 2, only 
models #15 and 16 are marginally acceptable. We can use these two models to set an upper 
limit to the amount of heavy element accretion during the early main sequence evolution of 
the Sun. 



6. Discussion 

There has been interest during the past few years in the effects of accretion on the Sun 
during its evolution. With this application in mind, the main purpose of this paper was to 
study the properties of a number of non-standard solar models in which the convection zone 
is richer in heavy elements than the interior, and to explore the constraints of helioseismology 
on such models. 
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A set of standard solar models were first constructed using the same physics and stellar 
evolution code as GD97, to serve as benchmarks for our non-standard models. Although our 
purpose was not to achieve the best possible fit to observation, we find that our standard 
solar model ^^20 satisfies the most stringent constraints from helioseismology at least as well 
as the best published standard solar model. It is just outside the la limit for the (Z/X)surf 
constraint, but agrees very well with the SNO Collaboration neutrino flux. It is apparent 
that a standard solar model in full agreement with all up-to-date observational data would be 
achieved with the same physics as in GD97, for a Zinit intermediate between model #19 and 
model #20, i.e. for 0.020<Zinit <0.022. It would also be interesting to probe intermediate 
values of Zj^t, with 0.80<Zint <1.00. 

For this work, we have only considered Z accretion from cometary material after the 
pre-MS phase, but what about the possibility of helium accretion as well, and increased 
accretion during an earher, pre-MS era? If most of the accretion takes place in the pre-MS 
phase of solar evolution, we must consider two phases. In the early phase, the proto-Sun is 
fully or nearly fully convective (the Hayashi phase). Since the accreted material is mixed 
efficiently in the convective region, composition gradients in the deep interior would not be 
expected. In the later phase of pre-MS evolution (the Henyey phase), the radiative core gets 
progressively larger until the main sequence is reached. This phase of evolution, which takes 
place on a thermal timescale, is relatively short compared to diffusion timescales, and could 
leave behind a composition gradient in the interior, below the convection zone. Since we 
do not understand the region just below the present convection zone well, we cannot rule 
out the possibility of some composition in-homogeneity in the outer radiative envelope due 
to late pre-MS accretion. However, two factors argue against the existence of significant 
composition gradients in the outer envelope: (1) the presence of shear-induced turbulence 
in the tachocline region, and (2) the close agreement of the observed sound-speed and the 
SSM sound-speed in the radiative region below the tachocline layer. 

None of our non-standard solar models agree as well with observation as our best stan- 
dard models. But two of the least extreme non-standard solar models, models #15 and 16, 
come close to satisfying all of the observational constraints. For this reason, they provide a 
realistic upper limit of ^2M^ of iron (or ^40M^ of meteoric material) to the accretion expe- 
rienced by the Sun during its early main sequence phase. This conclusion is compatible with 
Henney & Ulrich's (1998) earher null result that the accretion of SM® of meteoric material 
on the Sun could not at this point be detected by seismology because of the uncertainties. 
Similarly, a solar enrichment of O.SMq in iron, as suggested recently by Murray et al. (2001), 
is beyond detectability by seismic means at the present time. 

We are indebted to Prof. J. N. Bahcall for his advice and help with the neutrino 
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Table 1. Solar Model Characteristics 
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Table 2. Solar Model Nuclear Data^ 
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=■1-16 Non-Standard Solar Models, 17-20 Standard Solar Models 

''Values given in units of cm~^s^^, except <I>("^^C1) and <J>(^^Ga), which are given in SNU. 
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Table 3. Solar Model Constraint Values' 
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Fig. 1. — Grid of p-mode frequency difference plots, (model - Sun), for all models computed 
in this work. Lines connect data with a common /-value. Dashed lines correspond to 1=0-4, 
10 and 20. Sohd hues correspond to /=30, 40, 50, 60, 70, 80, 90, and 100. Each plot is 
annotated with the corresponding model number (see Tables 1-3), the surface Z/X ratio, 
and the neutrino flux, in cm~^s~^. Observational data used are from the MDI-SOHO 
data set (Rhodes et al. 1997). 
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Fig. 2. — p-mode frequency difference plot, (model - Sun), for model #20 computed in this 
work. Lines connect data with a common Z-value. Dashed hues correspond to l—O-A, 10 

and 20. Solid lines correspond to /=30, 40, 50, 60, 70, 80, 90, and 100. Listed are the 
surface Z/X ratio and the neutrino flux, in cm~^s^^. Observational data used are from 
the MDl-SOHO data set (Rhodes et al. 1997). The error bars indicate lOo" errors in the 
data, averaged over 250//Hz frequency bins, for all /-values included in this work. 
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Fig. 3. — Grid of small spacing difference plots, (model - Sun), for all models computed in 
this work. Lines connect data from a particular model, indicated by Zjnt and Zjnit. Each row 
contains models of the same Zjnt, where Zint=0.30Zinit includes models #1-4, Zint=0.50Zinit 
includes models #5-8, Ziiit=0.65Zinit includes models #9-12, Zint=0.80Zinit includes models 
#13-16, and Zint=Zinit includes the standard solar models #17-20. Each column contains 
results for a common Z- value, with ^=0 in column 1, /=1 in column 2, and /=2 in column 3. 
Observational data used are from the MDI-SOHO data set (Rhodes et al. 1997). 
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Fig. 4. — Grid of fractional difference plots for the sound-speed ((Sun - model)/model), 
for all models computed in this work. Solar sound-speed data were determined by Basu, 
Pinsonneault & Bahcall (2000) using the MDI-SOHO solar frequency data set (Rhodes et 
al. 1997). 



- 21 - 



Q. 





_ 1 1 1 1 1 1 1 


1 1 1 1 .1 1 1 1 1 






1 1 1 1 1 1 1 1 1 




.1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 






1 1 1 1 1 1 1 1 1 




0.2 


















































- 




























0.1 
























- 







U 1 






#<■ 




.III I'i'l 1 












0.2 


.III \ ] 1 


1 1 M 1 1 1 1 M 


-+ 


.III 1 ]"l 1 


1 1 1 1 1 1 1 1 1 1 


-H 




1 1 1 1 1 1 1 1 1 1 


-+ 


.1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 




























_ 


0.1 





























































HP, 






J/7 








- 


0.2 


.1 1 II 1 II 


1 1 1 1 1 1 1 1 1 1 1 


-+ 


.111111 


1 1 1 1 1 1 1 1 1 1 


-H 


.111111 


1 1 1 1 1 1 1 1 1 1 


-+ 


.111111 


1 1 1 1 1 1 1 1 1 




0.1 
























- 

















































M 1 2 




0.2 


.1 1 M 1 M 


1 1 1 1! M 1 1 1! 


-+ 


.1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


-h 


.1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


-+ 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


-+■_ 


0.1 



















































































#13 






#14 






#15 






#16 




0.2 
0.1 


.1 1 1 1 1 1 1 


11 1 11 1 1 1 1 11 


-+ 


.1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


-+ 


_l 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


-+ 


J 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


-+ 























































■, , . , 1 . . 


#17 
1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 


#18 

, 1 . , , . 1 , . 




1 1 1 1 1 1 1 


#19 
1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 


#20 
1 1 1 1 1 1 1 1 1 





0.25 0.5 0.75 0.25 0.5 0.75 0.25 0.5 0.75 0.25 0.5 0.75 

r/Ro 



Fig. 5. — Grid of fractional difference plots for density ((Sun - model) /model), for all models 
computed in this work. Solar density data were determined by Basu, Pinsonneault & Bahcall 
(2000) using the MDI-SOHO solar frequency data set (Rhodes et al. 1997). 
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Fig. 6. — Small spacing difference plot (model - Sun), for standard solar model #20 with /=0. 
Observational data used are from the MDI-SOHO data set (Rhodes ct al. 1997). The error 
bars plotted represent errors in the data. The observational uncertainties are approximately 
±0.07 /jMz for the small spacing determination. 
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Fig. 7. — Fractional difference plot for the sound-speed ((Sun - model)/model), for standard 
solar model ^^20 (zoom in of Figure 4). Vertical error bars indicate la errors in the inversion 
results due to errors in the data. Horizontal error bars are a measure of the resolution of the 
inversion. 




Fig. 8. — Fractional difference plot for density ((Sun - model) /model), for standard solar 
model ^20 (zoom in of Figure 5). 



